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Abstract. - We consider a mean- field system described by a general collective variable M, 
driven out of equilibrium by the manipulation of a parameter ^. Given a general dynamics 
compatible with its equilibrium distribution, we derive the evolution equation for the joint prob- 
ability distribution function of M and the work W done on the system. We solve this equation 
by path integrals. We show that the Jarzynski equality holds identically for these dynamics, 
both at the path integral level and for the classical paths which dominate the expression in the 
thermodynamic limit. We discuss some implications of our results. 



In a class of micromanipulation experiments a system initially in thermal equilibrium is 
subject to the external manipulation of a given parameter fi, and the work W done on it is 
evaluated [1,2]. This allows one to gather some interesting information on its structure, and 
in particular makes it possible to evaluate the free energy difference between the initial and 
final values of the manipulated parameter, by means of the Jarzynski equality [3] 

(e-^'^) = e-'^^^. (1) 

In this expression, (3 = l/ksT, where T is the temperature, the average on the Ihs is per- 
formed over the realizations of the process, and AF is the difference between the values of the 
equilibrium free energy F for the final and the initial values of the manipulated parameter ji. 

In order to evaluate correctly the Ihs of eq. (|l|) it is important to gather information on 
the probability distribution function of W for a given protocol In a previous work [4] we 
pointed out that a convenient way of achieving this goal is by considering the joint probability 
distribution function ^cr{W) of W and the microscopic state a of the system. We considered 
a comparatively simple system, like a pulled biomolecule subject to kinetic barriers, and 
obtained the corresponding work distributions. It is clear, however, that one should consider 
in general more complicated systems, described by collective coordinates and by more complex 
stochastic equations of motions. 

In a recent work [5] F. Ritort introduced a path integral representation of the work distri- 
bution for a very simple system, namely a collection of noninteracting spins, and showed how 
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its properties could be derived by drawing on thermodynamical analogies. Here we wish to 
generalize both approaches, by considering a system described by a general collective coordi- 
nate M, an equilibrium free energy J-^(M), depending on a parameter /i which is manipulated 
according to a protocol fi{t), < t < t{, and a general stochastic evolution described by a 
partial differential equation of the form 

dP 

-^-np. (2) 

where P{M, t) is the probability distribution function of the collective variable M, and H is 
a differential operator, depending on /i, and compatible with an equilibrium distribution 

Peq(M;^)oce-''^-(^'^). (3) 

We should have, therefore, for a given fi and any value of M, 

7^Peq(M;^) = 0. (4) 

During the manipulation of fi, the evolution of the system is described by a stochastic trajec- 
tory M{t), and the work done on the system reads 

w^f\,m'-^^^ (5) 

Jo d^i 

The evolution equation for the joint probability distribution $(M, W,t) of M and W can be 
derived by following the steps of ref. [4] : 

A similar equation for the joint probability distribution of the microscopic state and the work 
has been derived in ref. [6]. It is actually more convenient to introduce the corresponding 
generating function \I>(M, A,t) for the work distribution: 

■^{M,X,t) = j dWe-^^ ^{M,W,t). (7) 

One can readily see that ^ satisfies the following evolution equation: 

|_.«*-A,^*, (8) 

Thus, if A = /3, 

V{M, t) = exp (Af )] /Zo, (9) 

(where Zq is the partition function of the initial state) solves this equation and satisfies its 
initial condition 

«'(M, A = /3,0) = 7'(M,0) = — . (10) 

Thus '^{M,\=(3,t) = P(M, t), [7] and by integrating this function upon M one obtains the 
Jarzynski equality (||). We have thus shown that the solution of eq. (||) satisfies the Jarzynski 
equality identically. 
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In order to obtain one can attempt to solve directly eq. (|^) by numerical means. This 
approach runs into numerical difficulties when the system size N becomes large, because the 
corresponding distributions become narrower and narrower. Much insight is gained by looking 
at its path integral solution in the thermodynamic limit, as pointed out in ref. [5]. By going 
through the steps explained, e.g., in [8], one can see that, if the differential operator H is given 

by 

=E9M^{5.(M).}, (11) 

k=0 

one has 



n i.M{t) = M f n ^1 

«'(M,A,i)= / dMo / VjVM eiipl I At' C{t')\ *(Mo,A,0), (12) 

J Jm(0)=Mo iJ } 



f.M{t)=M 
'M(0)=A/o 

where the Lagrangian JZ is given by 



(^7M + H(7,M)-AA^) 



(13) 



'y{t),M{t),f,{t) 

Here 7i(7, M) is given in terms of the differential operator (O), 



W(7,M) = ^7'.9fc(M,Mi)), (14) 



fc=0 

and the path integral is defined by 

AT 



f^^.r 1- -fV d7(fcAt) dM(fcAi) 

AAAt=t k=l 

It can be readily seen that this formulation encompasses that discussed in ref. [5] . 

We denote by N the size of the system, and set M = Nm, J^^,{M) = Nf^{m), H{j,M) = 
NH{'j,m). As iV ^ oo, the path integral is dominated by the classical paths (7c(f), mc(t)), 
satisfying the equations of motion 

= m+— ; (16) 



dH f 
^ ~^ dm ^ dmdiJL ^ ^ 



The condition at i = for these equations of motion must now be discussed. For all values of 
A, ^'(Af, A, t—Q) satisfies eq. (|l0|). Thus the integral over the variables at t = takes the form 

*(M, A, i = 0) = = / ^° ° q1o(Mo-m,) ^ (^g) 

Zo J 2tt\ Zq 

Evaluating this integral by the saddle point method one obtains TOq = mo and that 70 and 
mo satisfy 

Q^r'7+|^ = 0. (19) 

Thus, for any given value of mp one can evaluate the corresponding value of 70 by imposing 
eq. (|l^) at t = 0, and then integrate the equations of motion (|l^,|l^) for < t < tf. This 
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will yield the values (rnf,7f) of m{t{) and 7(if) respectively. By evaluating the corresponding 
action one can thus obtain the value of ^'(TVmf, A, if). The generating functional of the work 
distribution r(A,tf) is given by the integral of '^{Nm{, X,tf) with respect to mf, see eq. (^. 
As pointed out in [5], this corresponds to evaluating ^'(iVmf , A, if) for that particular solution 
of the equations of motion which satisfies 7f = 0. Thus one must solve the equations of motion 
( P^ , p7| ) by imposing the boundary condition jig] ) at i = and the other condition j{t) = 
at i = if. This can be done by a shooting method: one solves the equations as a function 
of the initial value mo of m and chooses that value niQ for which 7(if) = 0. (This had not 
been necessary in ref. [5], since in that case the equation for 7 did not involve to.) This can 
be done more easily if one observes that, for A = 0, the solution j{t) = is the correct one. 
One can build on it to evaluate successively the solution for slowly increasing values of |A|. 
Alternatively, one may resort to a relaxation method, which is more convenient for slower 
manipulation protocols. Once the generating function r(A,it) of the work distribution is 
obtained, the distribution itself can be obtained by a saddle point integration, which amounts 
to a Legendre transformation on logF: 

P{W,tf) = J dXe^^+'°sm^)] = eX'w+ios[rix%u)] ^ (20) 

where A* is implicitly defined by ^Alogrj^, — —W. It is worth noting that we obtain in 
this way a probability distribution of the work per spin, w — W/N, which has the form 
P{w) oc exp [Ng{w)], where g{w) is a function of w alone. 

We now show that the solution obtained in this way satisfies the Jarzynski equality iden- 
tically. We first show that, for A = /3, the solution of the classical equations of motion (nq,0) 



satisfy an equation analogous to (19) at all times. By multiplying both sides of eq. by 
e"'''*^ and integrating by parts over M one obtains 

J dM H(7, M) e-/'^M(*^)-^^'^ = 0, (21) 

where 'H{j,M) is given by (|l4|). Evaluating this integral by the saddle point method in the 
large N limit, we obtain 

iJ(7,TO*)=0, (22) 



if 7 and m* are related by (|19|). By differentiating eq. ( P2D with respect to 7 at fixed /j, we 
obtain 

dm* 



dH dH 
dm 



= 0. (23) 



Let us also take the derivative of eq. ( [l9| ) with respect to 7 at fixed /u, we obtain 

By multiplying both sides of eq. (|23| ) by j^jdw? and substituting eq. (|2^), we obtain the 
following relation 

drn'- (?7 ^ dm ' ^ ^ 

which holds when 7 and m are related by eq. ( pg|). We can now evaluate the time derivative 
of the Ihs of eq. (p^), when 7 and m satisfy eqs. ( |l6|Jl7| ). We have 

Q ^ _ _ ^^f, = _ + - ^^/i . (26) 

dm^ dmdjjL \dm drnd/i J dm^ d^ dmd^ 
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The second and the last term cancel out. Substituting eq. (|25|), we see that also the first and 
the third therm cancel out. Thus if 7 and m satisfy eqs. ( |16| , |17D at all times, and satisfy eq. 
( |l9| ) at a given time, they satisfy this last equation at any time. 

Thus, for A = /3, the Lagrangian, evaluated along the classical path, is given by 



7m - f3^^— 



= N 



am Ofi 



(27) 



where we have exploited eq. (|19|). Substituting this expression in eq. (12) one recovers eq. 
and the Jarzynski equality. 

In the following we consider a simple mean-field Ising-like system, with free energy 

J 



^^^^^ " '^N^'^^ ~ ~ ^^^^'^^^ 



where 



S{M) = NkB 



In 2 



1 



, m 

ln(l — m") In , 

2 ^ ^ 2 V 1 



1 



(28) 



(29) 



The evolution of the system is governed by a Fokker-Planck equation with differential operator 



n- = 



d 
dM 



UJQ 



[dM 



(30) 



We show in fig. |l] the results of our calculation for this system, with (3 — 1, J — 0.5, 
N — 50, and where the external field h is manipulated according to the protocol 



Wo 



h{t) = ho + {hi-ho) — ; 0<t<tf. 



(31) 



We have set ho = —hi = —1, tf — 2. We compare the predictions with the results of the 
simulation of the Langevin equation underlying eq. (pO[), integrated according to the Heun 
algorithm [9]. We see that in this case the agreement between the simulation results and 
the distribution obtained by our analysis is quite good: nevertheless, one can also conclude 
that even for such a small system the possibility of using the results of the nonequilibrium 
manipulation protocol to obtain information on the equilibrium free energy is rather remote. 
Indeed, we have plotted in fig.|l| the quantity Q{w) = P{w)e~^^'^ , whose integral should 
give e~^^^ (which is equal to 1 in our case). We can see that not a single point among the 
simulated 10000 falls in the range in which f2(w) is essentially different from 0. 

When h = and (3 J > 1, our system undergoes a symmetry-breaking phase transition. 
Now assume that h — 0, (3J{0) < I, /3J{t[) > 1. Since jjU) = 0, the identically vanishing 
solution (7c = 0, TOc = 0) satisfies the equations of motion (|^ 
This would lead to 14^ = in contradiction with the Jarzynski equality. 



17) and the initial condition. 

However, this solu- 
tion does not correspond to the extremum of the exponent in eq. (^ij), which is given by the 
symmetry-breaking solution. It is therefore necessary to perform a quasi-average, by intro- 
ducing a small symmetry-breaking field h, evaluating the classical action, and going to the 
limit — > at the end of the calculation. We show in fig. || the results of this calculation. We 
have set ho — hi = 0.01, and manipulated J according to the protocol 



J{t) = Jo + {Ji- Jo)^, 



<t<ti, 



(32) 



with Jo = 0.5, Ji = 1.5, tf — 10, P = 1. Notice that for such a protocol the work done on the 
system is non-positive, since w — — J dtJ{t)m'^{t). The figure contains the results of 10000 
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Fig . 1 - Results for the system described by the diflerential operator ( po| ) with equihbrium free energy 
(p^, manipulated according to the protocol Continuous line: probability density P{w) of the 

work "per spin" w — W/N, with N=50. Dashed line: histogram of the work obtained by 10000 
simulations of the process via the Heun algorithm [9]. Dotted line: = P{w) exp{—PNw), whose 

integral verifies the Jarzynski equality. One can see that there are no simulation points in the range 
in which is essentially difi'erent from zero. Thus, in this case, a quantitative verification of the 

Jarzynski equality would be practically impossible. 



simulations each with = 1 (continuous line), = 10 (dashed line), iV = 100 (dotted line), 
and the expectations of the theory for N = 100 (dash-dotted line). The agreement between 
theory and simulations is poor for the smaller values of N, reasonable for N — 100, and 
improves as N grows (other data not shown). However, notice that for small values of N the 
simulated values of w are consistently smaller than those obtained asymptotically for large N, 
and lie closer to the thermodynamical value t«icv denoted by a vertical line. Indeed, for small 
values of h and large values of N the mean-field system has difficulties in keeping close to the 
true thermodynamic minimum, whereas a smaller system can reach it by activated processes. 

We have derived the differential equation (^ for the joint probability distribution of the 
collective coordinate M and the work W for a system evolving according to a general dynamics. 
We have shown that the solution of this equation satisfies the Jarzynski equality identically. 
We find that this equation can be solved by the path integral (p^). This solution can be 
evaluated, asymptotically for large system sizes N, by a saddle-point integration, leading to a 
two-point boundary value problem for the solution of the classical equations of motions, thus 
generalizing the approach of ref. [5] to interacting systems. We have seen that the Jarzynski 
equality can be effectively exploited to obtain the free energy difference only for very small 
systems. We have pointed out that the method has to be handled carefully if the manipulation 
protocol leads in the vicinity of a phase transition. (The situation is even more delicate if 
the protocol crosses a discontinuous transition line.) It is possible to extend the approach to 
disordered mean-field systems like the Hopficld model. By allowing one to explore the work 
probability distribution function, the formalism can be of use in identifying the manipulation 
protocol most likely to yield the desired results. 
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Fig . 2 - Results for the system described by the differential operator ( po| ) with equilibrium free energy 
(p8[), manipulated according to the protocol (|3^. Histograms over 10000 simulations for each value 
of iV. Continuous line: iV = 1. Dashed line: = 10. Dotted line iV = 100. Dash-dotted line: 
Calculated P{w) for N — 100. The calculated P{w) (not shown) for smaller values of N have the 
same general behavior as for N = 100, but are wider. Vertical line: Thermodynamic value of the 
work Wicv = AF/N. 

We are grateful to Felix Ritort for a critical reading of the manuscript. LP is grateful to 
Silvio Franz for thought-provoking remarks. 
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